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ABSTRACT 

Nuclear stellar discs (NSDs) can help to constrain the assembly history of their host 
galaxies, as long as we can assume them to be fragile structures that are disrupted 
during merger events. In this work we investigate the fragility of NSDs by means of N- 
body simulations reproducing the last phases of a galaxy encounter, when the nuclear 
regions of the two galaxies merge. For this, we exposed a NSD set in the gravitational 
potential of the bulge and supermassive black hole of a primary galaxy to the impact 
of the supermassive black hole from a secondary galaxy. We explored merger events 
of different mass ratios, from major mergers with a 1:1 mass ratio to intermediate 
and minor interactions with 1:5 and 1:10 ratios, while considering various impact 
geometries. We analyse the end results of such mergers from different viewing angles 
and looked for possible photometric and kinematic signatures of the presence of a disc 
in the remnant surface density and velocity maps, while adopting detection limits from 
real observations. Our simulations show that indeed NSDs are fragile against major 
mergers, which leave little trace of NSDs both in images and velocity maps, while 
signatures of a disc can be found in the majority of the intermediate to minor-merger 
remnants and in particular when looking at their kinematics. These results show that 
NSDs could allow to distinguish between these two modes of galaxy assembly, which 
may indeed pertain to different kinds of galaxies or galactic environments. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution, galaxies: 
nuclei, galaxies: interactions, galaxies: kinematics and dynamics 


1 INTRODUCTION 


Nuclear stellar discs (NSDs) were first discovered in im- 
ages taken with the Hubble Space Telesco pe (HST, e.g., 
Ijaffe et ^ 1 19941 : | van den Bosch et al.l 1 19941 ) , which allows 
for the detection of such structures that are often of only 
tens to hundreds of parsecs across an d which reside mainl y 
in the centre of early-type galaxies llPizzella et al.l [20021) . 
Previous studies have focused on the properties of few 
of such discs (e.g., MoreUi^taL I |2004 IPizzella et aLll200l 
IScorza fc van den Boschll 19981 ) but now thanks to the census 
of Ledo et al. I (120101 . hereafter Paper I) we have a sufficiently 
large sample to study these features more systematically. 

NSDs are not only interesting tools for determin- 
ing the masses of supermassi ve black holes (SMBH, 
Ivan den Bosch fc de Zeeuwlll99B ). but can also be a pow- 
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erful tool to trace galaxy merging history. In Paper I we 
have shown evidence pointing to the destruction of these 
discs during a merger event, even though if there is gas in¬ 
volved during such an event, part of this gas can fall to 
the centre and form a new stellar disc. Therefore, when 
present, NSDs provide us with a look-back time to the last 
merger experienced by their host galaxies, a precise one if 
they were born at that moment, or at least a lower limit if 
they formed afterward. This look-back time can be used to 
test and improve galaxy formation s cenarios and the predic¬ 
tions of semi-analytic al models (e.g.. lKhochfar fc SilkfioOGal : 
iDe Lucia et al.l |2006l ). In fact, NSDs are particularly apt 
for this role since their stellar age can be determined very 
precisely. Typical degeneracies between the stellar age and 
metallicity that affect population studies can indeed be re¬ 
duced using integral-field spectroscopic data and by know¬ 
ing in a dvance (thanks to a pho tometric disc-bulge decom¬ 
position [Scorz^)&j^end^[l9^) the disc-light contribution 
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(Sarzi et al., in preparation). What is still unknown to us 
is understanding more precisely to which extent these discs 
are fragile and what kind of merging events they actually 
trace. Do they trace only major mergers (as shown in the 
preliminary simulations published in Paper I), or do smaller 
satellite galaxies also have the ability to destroy them? If so, 
how small and in which circumstances? 

The impact of mergers has previously been explored 
in the case of kinematically decoupled cores (KDC, se e 
iBois et ^1201 if l or metallicity gradients (lKobavashill2004h . 
but so far little is known about the survival of NSDs. In 
fact, in this study we are primarily interested in establishing 
whether the presence of an NSD would still be detectable in 
the merger remnant by looking for the typical photometric 
and kinematic signatures of NSDs, such as discy isophotes 
and regular disc-like rotation, as observers would do on real 
galaxy images or velocity helds. At the same time, we will 
resort to our view of the disc particles at the end of our sim¬ 
ulations to better understand the disc disrupting process as 
a function of the parameters that we will explore, checking 
also for false positive detections. 

This work is organised as follows. In Sec. [2] we de¬ 
scribe the set up of our initial conditions and the various 
kinds of merging simulations that we considered. In Sec. [3] 
we describe the construction of surface density and velocity 
maps for our merger remnants as seen from different view- 
ing angles and our use of the technique called kinemetry 
llKrainovic et al.ll2006f l to analyse these maps and look for 
both photometric and kinematic signature of discs. In Sec. [4] 
we then assess how often such signatures would indicate the 
presence of a disc, and finally discuss our Hndings in Sec. [S] 


2 SIMULATION THOUGHT AND DESIGN 

We are interested in understanding how a galaxy merger 
would affect the innermost parts of a galaxy with a small 
nuclear disc embedded at its center. This is likely to depend 
on the nature of the encounter, for instance on it being either 
a major or minor merger, on the geometry of the collision, 
such as prograde as opposed to retrograde, and Hnally on 
the mass of the central disc, given that more massive discs 
should be more stable. In this paper we deal only with the 
magnitude and geometry of the perturbing event, and fo¬ 
cus on nuclear stellar discs, which can extend up to a few 
100 pc (Paper I) and can have stel lar masses comparable to 
that of their host galaxy SMBH JScorza fc van den Bosclil 
Il998h . Furthermore, we conducted a series of simulations 
meant to look only at the Hnal stages of a merging event, 
when the SMBH hole of the perturbing galaxy reach the 
central regions of the galaxy hosting the nuclear disc and a 
primary SMBH. Dynamical friction will indeed cause such 
a secondary black hole to sink towards the central regions 
of the primary galaxy or what would eventually become the 
merger remnant. Most of the stellar cusp around the sec¬ 
ondary is likely to be stripped during the early stages of 
the galaxy pairing, even though the total stellar mass that 
would remain bound to the secondary could depend on a 
plethora of factors, such as the initial mass concentration of 
the secondary galaxy or the mass ratio between the merg¬ 
ing galaxies, the orbital parameters of the mergers, an d so 
on I Callegari et al.l l2009l : IVan Wassenhove et al’l l2014[ l . At 


1:5 merger 



time (Myr) 


Figure 1. Time evolution for the separation between the primary 
and secondary black holes for the 1:5 mass ratio encounters. Dot¬ 
ted lines denote prograde encounters, whereas dashed lines repre¬ 
sent retrograde ones. The case of 1:1 (1:10) mergers would follow 
similar decay patterns although at shorter (longer) time scales. 


present we therefore decided to start by considering the sim¬ 
plest case of a naked secondary black hole. 

After setting stable initial conditions for a 10® Mq nu¬ 
clear disc with a radius of 100 pc that is embedded in a bulge 
with its own 1O®M0 SMBH, we introduce in the simulations 
an additional massive particle at an initial distance of 80 pc 
from the centre, on a circular orbit. This particle mimics the 
presence of a perturbing SMBH (the secondary, hereafter) as 
it is reaching the centre of the primary galaxy during the last 
stages of a galaxy merger. Inevitably, this secondary SMBH 
interacts with the star particles of the primary nucleus and 
sinks toward the centre while perturbing the distribution of 
the stars. After the secondary had reached the centre, the 
simulations were let to evolve for several dynamical times as 
to allow the system to relax to a point where we would not 
expect further perturbations to the dynamics of the stars. It 
is at this point that we will proceed to analyse the merger 
remnant, looking in particular to detect the same photomet¬ 
ric and kinematics signatures that observers would regard as 
indicative of the presence of a nuclear stellar disc. 

The simu lations were conducted using GADGET-2 
llSDringelll2nn!Tll . a parallel code designed for cosmological N- 
body/smoothed particle hydrodynamics simulations. Each 
simulation involved 600K stellar particles and two additional 
collisionless particles to model two SMBHs. The spatial reso¬ 
lution of our simulations is determined by the gravitational 
softening of the particles, which is equal to 0.5 pc for all 
the particles. 300K particles have been used to model the 
nuclear stellar disc of total mass Mdisc = 10® Mq and an 
initial outer radi us of 100 pc. The disc surface density, E, 
follows a Mestel llMestell[l963h proHle where E oc 1/R, and 
has an aspect ratio H/R « 0.05, similarly to previously used 
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BH 

Collision 

BH softening 

Mass ratio 

angle (°) 

(pc) 

1:1 

0 

0.5 

1:1 

30 

0.5 

1:1 

60 

0.5 

1:1 

90 

0.5 

1:1 

120 

0.5 

1:1 

150 

0.5 

1:1 

180 

0.5 

1:5 

0 

0.5 

1:5 

30 

0.5 

1:5 

60 

0.5 

1:5 

90 

0.5 

1:5 

120 

0.5 

1:5 

150 

0.5 

1:5 

180 

0.5 

1:10 

0 

0.5 

1:10 

30 

0.5 

1:10 

60 

0.5 

1:10 

90 

0.5 

1:10 

120 

0.5 

1:10 

150 

0.5 

1:10 

180 

0.5 


Table 1. Set of simulations with mass ratios, collision angles 
and the gravitational softening of the secondary SMBH. Collision 
angles 2 bh from 0 to 60° are prograde encounters whereas 2 bh 
values larger than 90° denote retrograde ones. 


relevant parameter space. This argument assumes a con- 
stant mass ratio bet ween the SMBHs and the host galaxies 
ijHaring fc Rixlliood h so that a 1:1 ratio in SMBH masses 
simulates a 1:1 galaxy merger. Even if this is the case at 
the beginning of a merger, the SMBHs can accrete a sig¬ 
nificant amount of gas, changing their initial mass ratio. 
Since the smaller secondary galaxy is more perturbed by 
the merger and suffers a higher gas inflow toward its centre, 
the secondary can accrete more (relatively to its in itial mass, 
ICalleeari et al.ir201ll : IWn Wassenhove et al.ll20l3 '). This re¬ 
sults in the formation of more equal mass SMBH pairs when 
they reach the nucleus of the merger remnant. For these rea¬ 
sons we consider our 1:10 encounters as providing a lower 
limit to the actual SMBH mass ratio in such small mergers. 
There is no consensus on what constitutes a major merger, 
which we have represented by a 1:1 mass ratio. Differ ent 
authors set limits as high as 1:2 (e.g., Broo k et al.| 2012|) or 
as low as 1:4 (e.g., iLop ez-Saniuan et ^ I I 2 OI 2 I : iBundv al.l 
I 2 OO 9 I : [Mailer et M]l2006l f . Therefore we chose to simulate a 
large minor merger, just below this limit, with a mass ratio of 
1:5 and a smalle r case of 1:10. For even more unequal merg¬ 
ers, works from ICalleeari et al. I (l201ll. I 2 OO 9 II suggest that 
such encounters do not lead to the formation of a close pair 
of SMBHs, since the tidal field of the primary quickly strips 
the secondary of most of its mass, thus decreasing the ef¬ 
ficiency of the dynamical friction at the early stages of the 
merger and leaving the secondary wandering at kpc separa¬ 
tions. 


nuclear discs (e.g., Dotti et ^l2007l : iDotti. Colni. fc Haardtl 
l200fil : lfecala et ^ 2005l i and t o the nuclear d iscs outcome of 
cosmological simulations fe.g.. iLevin 3 l2008ll . Another 300k 
particles have been used to model the central regions of the 
spherical stellar bulge hosting the disc. The bulge follows a 
Plummer profile 


3 Mbulge 
An 



( 1 ) 


where b is the core radius (6 = 50pc), r the radial co¬ 
ordinate and Mbulge the total mass of the spheroid, with 
Mbulge = 6.98Mdisc- Finally, the bulge and disc stellar struc¬ 
ture hosts a 1O®M0 SMBH. While the SMBH is initially 
at rest at the centre of the structure, the velocities of the 
star particles modelling the disc and the bulge are chosen to 
set the system in dynamical equilibrium. Since an analytic 
formulation of the distribution function of the composite 
structure is not available, we let the system to stabilise into 
dynamical equilibrium in isolation. We further verified that 
over the typical timescale of our simulated encounters (a 
few ten Myr, see Fig. [T]) the stability of such initial condi¬ 
tion was maintained. Individual simulations run time will 
vary slightly, depending on how fast the secondary BH sinks 
toward the centre and forms a binary with the primary (see 
Fig. ID). After this happens we let the simulations continue 
for a few (« 5) dynamical times to let the system relax. 
Only then we analyse the results. 

We designed and ran a set of simulations wide enough 
to probe a variety of merging scenarios with different mass 
ratios and orbital inclinations (see Tab. [T|. In these runs 
we vary the mass ratio between the primary and sec¬ 
ondary SMBH from 1:1 to 1:10 to explore the cosmologically 


3 SIMULATION ANALYSIS 

The simulations follow in time steps the positions and ve¬ 
locities of the black holes and of the disc and bulge stars. 
Once the system has relaxed and the distance between the 
primary and secondary black holes has stabilised itself we 
proceed to check whether the merger remnant still shows 
some photometric and kinematic signature of a nuclear disc. 
For this we will produce maps for the surface density and 
velocity field of the merger remnant, and analyse such maps 
using KINEMETRY (Krajnovic et al. 2006) in order to detect 
in particular discy isophotes and regular disc-like rotation. 
To help our interpretation of any photometric or kinematic 
detection or non-detection, we will also produce and visu¬ 
ally inspect maps for the surface density of the disc particles 
only at the end of our simulations. 


3.1 Stellar surface brightness and velocity maps 

The simulation outputs were first binned in regular grids 
to create the equivalent to surface brightness and velocity 
maps. This is a requirement of the code we used to detect 
the pre sence of discs, as it is of others such as IRAF task 
ellipse ll Jedrzeiewskil [l9^ 1 . Such grids also allow us to vi¬ 
sually follow the merging event and have an equivalent to 
an image as it would be the case in real observations. For 
simplicity, the grids were made only for the central regions 
where the disc was initially present, measuring 400 x 400 pc, 
using a grid spacing of 5 pc and a Gaussian spatial smooth¬ 
ing of the same size. We produced grids for both the entire 
remnant and for the disc particles only, even though the lat¬ 
ter maps will only serve for a qualitative assessment of the 
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Figure 2. Surface mass density (left panels, in Mq pc“^ and 
in logarithmic scale) and velocity fields (right panels in km s"*^) 
for the initial conditions of our simulations, showing, from top 
to bottom the disc particles, the bulge particles and the entire 
system. This figure zooms in a 400 x 400 pc region around the 
centre that encompasses the entire disc, which is being viewed at 
a 60° inclination. 


Figure 3. Kinemetry results for the initial conditions shown in 
FigO The blue, red and black lines show radial trends for the 
surface brightness, 04 , velocity and A:5/fel of the disc, bulge and 
all particles, respectively and as indicate at the top of each panel. 
A peak in the a 4 and low values of fcS/fcl indicate the presence 
of a disc. 


disc disruption. Because, as we will discuss later, the view¬ 
ing angle i affects our ability to detect a disc amidst the 
bulge light, we generated maps for cos(i) = 0.0, 0.1, 0.2, 0.3, 
0.4 and 0.5, spanning from an edge-on view to 60°. There is 
little or no point to a nalyse project i ons w ith values of cos(t) 
larger than 0.6 since iRix fc White! (Il990f) have shown that 
at this point it becomes very hard to photometrically detect 
the presence of discs. As a first example of these maps Fig. [2] 
shows the disc, bulge and total surface density maps with 
corresponding velocity maps, for the initial conditions of our 
simulations. 


3.2 Kinemetry 

A visual inspection of the maps for the surface brightness 
and velocity fields of the merging results allows already to 
appreciate the effects of the different merger simulations, but 
to better study the merger remnant and quantify whether 
a dis c would still be detect able in it we have used kineme¬ 
try llKrainovic et al.ll200(j ) to analyse the moments of the 
line-of-sight velocity distribution (LOSVD). The gridding 
described in the previous section was done not only to aid 


visualisation and qualitative assessment but also as a neces¬ 
sity for the photometric side of kinemetry. 

The KINEMETRY algorithm works under the assumption 
that the even moments of the LOSVD, such as the surface 
brightness, are described by a constant profile along ellipses, 
and that the odd moments, such as mean velocity, can be 
described by a cosine law along such contours. In partic¬ 
ular, we will be looking for peaks in the coefficient of the 
fourth cosine term on the Fourier series expansion of the de¬ 
viations from the elliptical isophotes (04) in the photometry, 
and for very small values of the k5/kl ratio (analogous to the 
04 ) in the velocity fields. Large 04 values will indicate discy 
isophotes, whereas small k5/kl values will ensure the pres¬ 
ence of a regular, disc-like rotation. As done in other pho¬ 
tome tric tools such as the IRAF task ellipse (Ijedrzeiewskil 
Il987l) . KINEMETRY begins by identifying the best-fitting el¬ 
lipses to the isophotes to retrieve the radial profiles for the 
average surface brightness, position angle (PA), flattening 
(g) and higher moments such as the 04 . Similarly, from the 
odd moments, we obtain a velocity curve (traced by the kl 
parameter) , a kinematic PA and flattening, and the k5/kl 
values (see lKrainovic et al.ll2006l . for more details). 

In Fig.[3]we show the main outputs from kinemetry for 
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Figure 4. Final total surface density maps and velocity fields (left panels) and correspoding kinemetric profiles for the 04 parameter, 
the velocity V and the k^jkl ratio (central and right panels) resulting from 1:1, 1:5 and 1:10 mass ratio encounters (top, middle and 
lower rows) with ibh = 30° and as viewed from an inclination i = 60°. In the panels showing our kinemetry results, the black lines 
show these measurements for the entire system whereas the blue lines are for the disc particles only. Vertical dashed lines indicate the 
original extent of the disc. We compute average values for those kinemetry parameters inside this radius of 100 pc, except for the most 
edge-on projections in the case of V and k5/kl, where we restricted our analysis to the inner 50 pc. 


the same initial conditions shown in Fig. [2] As expected, we 
can see 04 and velocity peaks in both the disc and total sur¬ 
face brightness (blue and black lines respectively) whereas 
the bulge (red line) stays flat. On the other hand, the k5/kl 
values for the disc and total velocity field present very low 
values, tracing disc-like rotation, whereas the bulge displays 
a noisy k5/kl profile. These are the quantities that will allow 
us to tell, once compared to appropriate detection thresh¬ 
olds (Sec. 3.3) if the photometric and kinematic signatures 
of a disc are present or not when looking at the final stages 
of the interaction. 


Fig. |4] shows examples of the surface brightness and 
velocity maps with corresponding kinemetry outputs for 
the final stages of the 1:1, 1:5 and 1:10 mass ratios after a 
iBH = 30° encounter and seen at a 60° inclination, similar to 
what was shown for the initial conditions. The bulge results 
have been omitted for simplicity as they are not key to our 
goals. Similarly, we do not follow the behaviour of the disc 
alone (blue lines) since, although interesting, what observa¬ 
tions trace is the surface brightness and average velocity of 
all the stars. For the purpose of detecting the signature of a 
disc, we will therefore look at and analyse the surface den¬ 
sity and kinematics of all the particles in our simulations 
(black lines), in particular, the region where the disc was 
originally present, which corresponds to the solid line in the 
plots and is delimited by the vertical dashed line, placed at 
the 100 pc radius. 


3.3 Disc kinemetric signatures 

In order to assess at the end of which kind of merger events 
we could still find the photometric or kinematic signatures 
of a disc, we considered the results of the kinemetry anal¬ 
ysis inside the central 100 pc region originally occupied by 
the disc, although for edge-on projections this radius was 
decreased to 50 pc when considering the velocity V and 
the k5/kl parameter, due to convergence problems with the 
KINEMETRY procedure. In Fig. [4] such a limit is indicated 
by a vertical dashed line, and the kinemetry descriptions 
inside this radius are plotted as solid lines. 

To identify the presence of a disc, we will use the mean 
values for the kinemetry parameters in such central re¬ 
gions, in particular to isolate fluctuations caused only by a 
fraction of the disc particles which had acquired extreme be¬ 
haviours, such as those ejected at large radii. Ideally, when 
computing such average values one would first guess the size 
of the disc by locating a clear peak of the 04 parameter, 
which indeed is a good gauge of the disc extent (see Fig. 5 
in Paper I). This approach would not be applicable to all our 
merger remnants, however, in particular for the most con- 
spicuos 1:1 mergers where the 04 signal is weak and noisy. 
On the other hand, when a disc signature is detected, we 
always found an 04 peak close to the original edge at 100 
pc of our initial disc, since in fact when the disc particles 
remain in a disc-like configuration, the extent of such disc 
remnants is similar to that of our initial conditions. This 
can be appreciated in Fig. [5] where we collected the surface 
density distribution of the disc particles at the end of all 
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Figure 5. Surface density for the disc particles at the end of our simulations. From left to right we show the disc remnants of 1:1, 1: 
and 1:10 mergers, both from an viewing angle of 60° and as seen edge-on, whereas from top to bottom we display simulations wit 
increasing impact 2 bh angles. The maps adopt the same logarithmic scaling and units as shown in Fig. [2] 
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our simulations as seen both edge-on and from the lowest 
inclination that we considered, of 60°. 

In the next section we will compare the average val¬ 
ues for the 04 parameter, the rotation velocity V and for 
the k5/kl ratio, with corresponding threshold values that 
indicate the presence of a disc. For the 04 parameter we 
considered a generous minimum average value of 0.02 for 
detecting discy isophotes in our simulations, which is com¬ 
parable to what is observed in real NSD-host ing galaxies 
(see e.g.. lMorelli et al.ll2004l : |Pizzella et al.ll2002h . As regards 
kinematics, we require an average velocity above lOkms”^ 
to detect global rotation, in line with typical errors on this 
quantity, and average k'bjkl values below 0.04 to recognise 
the presence of disc-like motions. The latter threshold, which 
lies just above the typical av erage value found in ou r ini¬ 
tial conditions, is taken from iKrainovic et al.l (l201l[) who 
based their judgement on integral-held da ta for 260 early- 
type galaxies from the ATLAS®^ survey (|CaDDellari et al.] 

Fig. [ 6 ] presents the central average values for the 04 , V 
and k5/kl parameters measured in the surface-density and 
velocity maps for the hnal step of each of our simulations as 
observed from different viewing angles from edge-on down to 
a 60° inclination. The different panels in these hgures group 
together simulations with similar impact angle Ibh. In each 
panel the size of the circles show the average 04 , V and 
k5/kl values, where the size of these symbols correspond 
to the mass-ratios of the merger events (i.e., with largest 
circles showing 1:1 mergers). Finally, the asterisks in Fig. [ 6 ] 
also show the average values of these parameters from our 
initial conditions. These are the figures that will be discussed 
in the next section. 


4 RESULTS 

Photometrically, the search for the presence of nuclear discs 
generally starts with the visual inspection of unsharped im¬ 
ages or structure maps, w hich highlight small scale features 
(e.g.. |Pizzella et ^l2002l . Paper I). This is then followed by a 
more quantitive analysis of the shape of the stellar isophotes, 
looking in particular for positive peaks in the 04 parameter. 
With integral-field spectroscopic data it will further be pos¬ 
sible to look for rotating structures that may not always 
show up photometrically, which are likely to be discs when 
Ending very low values for the feS/fcl kinemetric parameter 
that would indicate a regular, disc-like kinematics. 

4.1 Visual inspection of disc remnants 

Before analysing our kinemetric results on the global out¬ 
come of our simulations, we produce maps for the surface 
density of the disc particles only and inspect such disc rem¬ 
nants. This we do in order not only to visually assess the 
disruption of the disc but most important to gain some in¬ 
sights on our kinemetric Endings and the detection of disc 
signatures. As one would expect. Fig. [S] shows that major 
mergers (with a 1:1 mass ratio, left) have the most dra¬ 
matic impact on the disc, as its particles get scattered at 
large radii and vertically leading to lower surface density 
values and thicker remnants compared to intermediate and 
minor mergers (with 1:5 and 1:10 mass ratio, center and 


right). Fig.[5]further indicates that prograde encounters (top 
panels) affect the disc structure more than their retrograde 
counterparts (lower panels), which can be noticed by look¬ 
ing in particular at the 1:5 merger disc remnants. This vi¬ 
sual inspection also reveals interesting structures in the disc 
remnants, such as circumnuclear rings or the accumulation 
of stars in central clusters. 

4.2 Incidence of disc kinemetric signatures 

Having deEned earlier in Sec. I3.3l the thresholds for the de¬ 
tection of a disc in our simulations based on average values 
for the 04 , V and fc5/fcl parameters, and after checking that 
all our simulations reached a post-interaction relaxed state, 
we now proceed to assess in which kind of simulation and 
down to which inclination we can still End clear photometric 
and/or kinematic signatures of a disc. 

Starting from the plots for the 04 parameter (top panels 
in Fig. O, we re-observe the trend whereby the photometric 
signature of a disc gets more easily erased during the interac¬ 
tions involving a more massive secondary black hole. Indeed, 
1:1 mergers leave little or no trace of a nuclear disc in 81% 
of the cases whereas the Enal steps of the 1:5 and 1:10 merg¬ 
ers display average 04 values above the threshold suggestive 
of a disc structure in 73% and 86 % of the cases. Prograde 
encounters would appear to be generally more destructive, 
as can be noticed by comparing the trends for the 1:5 and 
1:10 remnants for ibh = 30 and 60° with their counterparts 
at iBH values of 120° and 150°, or by observing that most of 
the 1:1 merger remnants showing discy isophotes are found 
in the case of a perfectly retrograde encounter (Ibh = 180). 
In fact, the edge-on projection of the 1:1 prograde merger 
(ibh = 0 ) owes its large 04 values only to a distinct ring-like 
structure in the disc remnant (Fig. [SJ, which would not be 
mistaken for a disc upon a simple visual inspection. 

The difference between prograde and retrograde can 
be understood considering that stars are more easily de- 
Eected by the secondary black hole when this travels along¬ 
side rather than against them, which is also why dynam¬ 
ical friction on the secondary is more efficient in the pro¬ 
grade cases making it sink faster towards the centre. This 
can be appreciated back in Fig. [T] which shows the evolu¬ 
tion of the black holes separation as a function of time and 
where the prograde tracks for ibh = 0, 30 and 60 indicate 
always a faster sinking tha n their retrograde coun terparts 
for isH = 180, 150 and 120. iMeiron fc Laoil (l2013l ~) provide 
another example of this behaviour in the somehow differ¬ 
ent context of a binary black-hole interaction within a non¬ 
rotating bulge, which ends up showing a central rotating 
structure as the coalescing black holes evacuated preferen¬ 
tially the stars that co-rotate with them. The larger impact 
of a prograde encounter on discs was also already shown 
by IVelazauez fc Whit3 lll999h in the case of Milky Way-like 
galactic discs and satellites. 

Moving on to the trends for the peak velocity V traced 
by our kinemetric Et (middle panels in Fig. | 6 }, we observe 
that - except for the planar and polar 1:1 encounter and the 
planar prograde 1:5 merger - some level of bulk stellar rota¬ 
tion would always remain detectable in the remnant, even in 
cases where no photometric signature of the disc was iden- 
tiEed in the 04 proEles. This suggests, as one would expect, 
that while perturbing the radial and vertical structure of a 
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Figure 6. Average values for the a 4 parameter (top row), velocity V (middle row) and the kbjkl ratio (lower row) as measured at the 
final step of each of our simulations and as a function of our viewing angle (cosi). The different panels from left to right show the results 
from simulations of similar values for the collision angle ibhi with the size of the circles in each panel tracing the average a^, V and 
kbjkl values for simulations with increasing merger mass ratio (i.e., from 1:10, 1:5 to 1:1). As a reference, in each panel the asterisks 
indicate the corresponding values as measured in our initial conditions, whereas the dashed horizontal line indicates our threshold for 
detecting the photometric or kinematic signatures of a disc in the merger remnants (i.e., discy isophotes, coherent rotation and disc-like 
kinematics in the case as indicated by the 04 , V above and k5/kl below the thresholds, respectively). 


thin disc is relatively easy, entirely erasing the signature of 
the rotation of its stars is much more difficult, and requires 
a rather dramatic scattering of the disc stars at larger radii. 
This is for instance the case of the aforementioned ring-like 
structure in the 1:1 planar prograde merger, which is all that 
is left of the disc and does not contain a sufficient number 
of particles to imprint a signature of bulk rotation in the 
merger remnant, even when seen edge-on. 

Not all rotating remnants show regular disc-like kine¬ 
matics, however. Looking finally at the average values for 
the A:5/fcl parameters (lower panels in Fig. [6]) shows in par¬ 
ticular that nearly half of the 1:1 remnants do not rotate reg¬ 
ularly when they do show bulk rotation. On the other hand, 
except for a few projections of the prograde encounter, the 
remnants of 1:5 interactions generally display regular kine¬ 
matics, with this being always the case of the 1:10 remnants. 

Overall, a sufficiently ample and regular rotation is 
found in 38%, 72% and all of the remnants of the 1:1, 1:5 
and 1:10 mergers, respectively. Many of these also show pho¬ 
tometric disc signature that would suggest the presence of a 
rather thin disc component, whereas conversely the absence 
of a photometric signature, which occurs in particular in 1:1 
mergers, could lead to infer the presence of a rather thick 
disc (as Fig. [6] indeed shows). Based on this limited number 
of simulations we can draw an initial general picture for the 
chance of nuclear stellar discs surviving a galaxy merger and 
of being detected, considering in particular that co-planar 
encounters are, for a purely probabilistic perspective, much 
less likely to happen than merger events where the secondary 
black hole initially moved in an highly inclined orbit. We can 
account for this effect by correcting for a sin(iBH) factor our 
raw fractions of remnants with kinematic and photometric 
disc signatures, which effectively removes the perfectly co- 



fphot 


fkin 

Mass ratio 

raw 

sin(iBH) 

raw 

sin(iBH) 

1:1 

19.0 

16.1 

38.1 

45.5 

1:5 

61.9 

72.8 

71.4 

80.1 

1:10 

85.7 

85.0 

100.0 

100.0 


Table 2. Fractions in percent of merger remnants where the disc 
can still be detected, for the three mass ratios and according to 
whether the disc is detected photometrically or kinematically. The 
“raw” columns list uncorrected fractions, whereas the “sin(iBH)” 
columns list fraction that are corrected by this factor to take into 
account that co-planar encounters are much less likely to happen 
than those with higher ibh inclinations. 

planar cases and accentuates the impact of polar or nearly 
polar impacts. Since the former simulations tend to be the 
more disruptive for the disc, this correction tends to produce 
larger detection rates for discy or regularly rotating struc¬ 
tures. These fractions do not vary dramatically, however, in 
particular since the perfectly polar impact can also be quite 
damaging for the disc. The values for both our raw and 
sin(iBH)-corrected fractions of merger remnants with photo¬ 
metric discy features or regular velocity fields are reported 
in Tab. [21 


5 DISCUSSION 

Using a suite of numerical simulations we have explored the 
fragility of nuclear stellar discs in order to establish their use¬ 
fulness as tools to constrain merging histories of early-type 
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galaxies, looking in particular at the same photometric and 
kinematic signatures that observers would regard as indica¬ 
tive of the presence of a NSD. Our results show that NSDs 
are destroyed in major merger situations whereas photomet¬ 
ric and, even more generally, kinematic signatures of a disc 
are still detectable in the remnants of minor mergers (see 
Tab.IH). Integral-field studies should therefore prove partic¬ 
ularly useful to compile a more complete census of NSDs 
than already done based solely on HST image s in Paper I. 
In pa rticular, new instruments such as MUSE llBacon et alJ 
l20ldl . with its superior sensitivity and the ability to work 
in conjunction with adaptive optics, will be very useful for 
detecting NSDs. Additionally, MUSE should also prove an 
excellent tool for dating the age of these structures, improv¬ 
ing on the first experiment of Sarzi et al. (in preparation) 
based on VIMOS data. 

A large sample of NSDs with well constrained values for 
their stellar age, combined with the knowledge gained here 
on the fragility of such structures, should allow to paint a 
more detailed picture for the assembly history of galaxies. 
For instance, finding evidence for a NSD both kinematically 
and photometrically - that is - finding the signature of a thin 
disc, would suggest according to our analysis that there is 
only a ~ 20% chance that the host galaxy of such a disc 
would have experienced a major merger since the formation 
epoch of the disc, as indicated by its stellar age. In fact, 
our estimated incidence of nuclear-disc hosting galaxies that 
would still show evidence for a disc after a merger are most 
likely conservative limits, in particular since we have not 
considered here the possible presence of a stellar cusp around 
the secondary incoming SMBH, and the additional negative 
impact that would come with such an added stellar mass. 

Galaxy evolution models p redict merging trees of a 
combination of mer gers (e.g., iKhochfar fc Sillj l2006al [bl: 
iGonzalez et al.l I 2 OIII I and there is a great debate on the 
influence of major and minor mergers in the size and mass 
growth of gala xies, with current theory favouring 1:5 - 1:10 
merg ers (e.g., lOser et ^ 1 I 2 OI 2 I : iBedorf fc Portegies Zwarj 
l2013h . This corresponds to the minor merger range of mass 
ratios studied in this work. Yet major mergers cann ot be 
discarded, as other studies (e.g.. Ilhdman et al.ll2013l l have 
suggested that in recent epochs they may have caused large 
mass increases. Therefore, the ability of these NSDs to dis¬ 
criminate between merger scenarios is central to today’s 
questions. 

We are aware that the set up of the simulations pre¬ 
sented here is simplified. For instance considering a Mestel 
profile for the nuclear disc and a purely spherical Plummer 
distribution for the bulge does not convey the complexity of 
real nuclei where e.g., the disc is better described by an ex¬ 
ponential law and bulges can be flattened and present both 
core and cuspy profiles. Yet, with a bulge mass approxi¬ 
mately 5 times that of the disc within 100 pc, the corre¬ 
sponding disc to total surface brightness ratio of our initial 
mode l does span the ob served values between 20 and 40% 
(e.g.. iMorelli et al]l2004h . depending on the inclination. We 
therefore regard this work a first exploratory study, which 
could be improved on the one hand by exploring a larger 
parameter space, for instance by varying the mass of the 
nuclear disc compared to its host bulge and central SMBH, 
or on the other hand by considering a more realistic model 
with a stellar cusp around the secondary black hole. The 


actual amount of stellar mass bound to the sinking nucleus 
of a satellite galaxies may well remain a free parameter in a 
future analysis, unless new insights will emerge from large r 
sets of simulations such as those of ICallegari et aP (l2009li . 
but in general we expect that the presence of a stellar cusp 
around the secondary will further contribute to the disc dis¬ 
ruption. 

To conclude, we note that better simulations may as 
well shed more light also on the origin of the kinematically 
decoupled cor es (KDG) that are often found in early-type 
galax i es (e.g.. Bender 198^ fianx. Illingworth, fc HeckmanI 


Il989l : iMcDermid et al.ll2006ll . Indeed a significant number 
of our simulated encounters result in remnants that display 
regular rotation, characterised even by a low average /c5/fcl 
value, despite not showing any photometric a 4 signatures of 
discy structures. This resembles the situation observed for 
th e KDCs of massive e lliptical galaxies (see, e.g.. Fig C4 
of iKrainovic et ^1201 ill , where indeed a distinct kinematic 
central pattern is never accompanied by a significant change 
in the isophotes orientation, flattening or shape in the same 
regions. To us, this similarity suggests a possible link be¬ 
tween KDCs and puffed-up thick discs, and a closer look 
to the higher-order moments of the stellar LOSVD, both 
in future simulations and in high-quality integral-field data, 
should provide more valuable insights on the nature of these 
structures. 
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